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Abstract: Polydimethylsiloxane (PDMS) is commonly used as the coated polymer in the 
solid phase microextraction (SPME) technique. In this study, the partition coefficients of 
organic compounds between SPME/PDMS and the aqueous solution were compiled from 
the literature sources. The correlation analysis for partition coefficients was conducted to 
interpret the effect of their physicochemical properties and descriptors on the partitioning 
process. The PDMS-water partition coefficients were significantly correlated to the 
polarizability of organic compounds (r = 0.977, p < 0.05). An empirical model, consisting 
of the polarizability, the molecular connectivity index, and an indicator variable, was 
developed to appropriately predict the partition coefficients of 61 organic compounds for 
the training set. The predictive ability of the empirical model was demonstrated by using it 
on a test set of 26 chemicals not included in the training set. The empirical model, applying 
the straightforward calculated molecular descriptors, for estimating the PDMS-water 
partition coefficient will contribute to the practical applications of the SPME technique. 

Keywords: polydimethylsiloxane (PDMS); solid phase microextraction (SPME); 
partition coefficient; polarizability; molecular connectivity index 
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1. Introduction 

Solid phase microextraction (SPME) is a solvent-free sample preparation technique. The outer layer 
of the fiised silica rod of the SPME device is coated with poljmieric materials, such as 
polydimethylsiloxane. When the SPME device is placed in a sample matrix, the analyte is extracted 
and sorbed from the sample matrix onto the polymeric coating stationary phase. The analj^e is then 
desorbed from the SPME device and the concentrated extract is analyzed using an instrument such as a 
gas chromatograph. The SPME technique simplifies the four steps of sampling, extraction, 
condensation and introduction of the sample into the analytical instrument [1]. In comparison to the 
traditional solid phase extraction methods, the advantages of SPME are lower cost, easy handling and 
shorter time procedures. Therefore, SPME has been widely used in the analysis of many organic 
compounds in water [2,3]. 

Sorption of the analyte between the sample matrix and the coating phase may be assumed to be an 
equilibrium partitioning process and may be modeled by a linear isotherm as follows: 



Kfs=f (1) 



where C/and Cs are the equilibrium concenfrations of the analyte in the coating and the sample mafrix, 
respectively (M/L ); and K/g is the coating-matrix partition coefficient. 

For an aqueous mafrix, the total mass of analj^e in the coating and the sample mafrix can be 
represented as: 

C,V,=C,V,+cMKfsVf+ K K (2) 

where Co is the initial concentrations of the analyte in the sample matrix (M/L^); and V/and Vs are the 
volumes of the coating and the sample matrix, respectively (L^). 

The mass of analyte sorbed by the coating of SPME, Mf(M), can be determined as: 



M,=C,VrK,C/r^^^ (3) 



Prior to using the SPME technique for sampling and analysis, the first step is to select the 
appropriate poljmier as the SPME coated material to sorb the analyte. As shown in Equation 3, the 
mass of analyte sorbed by SPME is determined by the coating-mafrix partition coefficient Kfs. 
The proper application of SPME is primarily dependent on the partition coefficient of analyte between 
the polymeric coating and sample matrix. Therefore, there is a need to estimate the partition coefficient 
Kfs for determining the sensitivity of SPME exfraction [4—6]. 

Factors affecting the sorption of organic compounds onto the solid phase include electrical 
attraction, chemical affinity of the organic molecule, van der Waal's force and the hydrophobic nature 
of the organic compound [7]. The physicochemical descriptors, such as polarizability, molecular 
connectivity index, and indicator variable, have been found to be successful in modeling the 
parameters for partition processes like octanol-water, octanol-air and air-water partition coefficients as 
well as adsorption coefficients and aqueous solubility of organic compounds [8-12]. 
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This study is an attempt to develop predictive equations for the partition coefficients of a variety of 
chemicals in aqueous solutions onto SPME/PDMS based on their physicochemical properties and 
molecular structures. This was achieved by using the PDMS-water partition coefficients, Kj\v, retrieved 
from the literature [13,14] and relating them to their physicochemical descriptors, such as 
polarizability (O), molecular connectivity index (^x)? indicator variable (7), water solubility {WS), 
molecular weight (Mw), and octanol-water partition coefficient {log Kq^). The correlation analysis of 
Kjw can provide an understanding in the partitioning process between SPME/PDMS and aqueous 
samples. Further, the predictive equations will provide a means of determination Kfs for other 
chemicals, and facilitate the application of the SPME technique. 

2. Results and Discussion 

2.1 Correlation Analysis 

PDMS is a hydrophobic and nonpolar polymer. Therefore, the nonpolar attractive force causing 
physical adsorption of organic compounds onto PDMS surfaces is the van der Waals force. One of the 
attractive potential energies for the van der Waals forces is the dipole-induced dipole interaction, also 
called London dispersion force, between PDMS and organic molecules. When SPME/PDMS is placed 
in the aqueous solution, the dipole moment of PDMS creates an electric tield which polarizes the 
charges on organic molecules. The magnitude of induced change in dipole moment is determined by 
the polarizability of organic molecules [15]. As presented in Table 1, Kfw was found to correlate best 
with the polarizability of organic compounds (r = 0.911, p < 0.05). Therefore, the polarizability can be 
a good basis to understand the partition between SPME/PDMS and organic compounds in the 
aqueous solution. 

Table 1. Correlation analysis of Kjw for the training set. 



Descriptor Correlation 





/og = 1 .933 + 1 .026 ^ (r = 0.977, < 0.05) 


log Ko„ 


log = 0.3 1 9 + 0.88 1 log ^ow {r = 0.11 A, p < 0.05) 


\ 


log Kf„ = 0.750 + 0.862 'i (r = 0.743,^3 < 0.05) 


I 


log Kj„ = 4. 144 - 0.846 1{r = 0.583,/7 < 0.05) 


WS 


log Kf„ = 3.936 - 0.002 WS (r = 0.559, p < 0.05) 


Mw 


logKj„ = 1.859 + 0.017 M\v (r = 0.474,/; < 0.05) 



In general, the van der Waals forces of molecules increase with an increase in their molecular 
weights [7]. Table 1 shows that K^v was proportional to the molecular weights of organic compounds 
but only with a correlation coefficient r = 0.474 (p < 0.05). If the organic compounds were grouped in 
different categories, i.e., alkanes and aromatic hydrocarbons, for correlation analysis. Figure 1 
indicates that Kfw was signiticantly correlated to their molecular weights (r > 0.916). This result may 
imply that difflision of organic molecules in the PDMS fiber was not a governing mechanism for the 
mass transfer between PDMS and organic compounds. 
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Figure 1. Correlations of for the training set with their molecular weights. 
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According to the principle of "like dissolves like", it is plausible that a non-polar compound is more 
easily sorbed by PDMS. In general, the larger the octanol-water partition coefficient of a chemical, the 
lower the polarity. Several researchers indicated that log Kj^, of SPME/PDMS linearly increased with 
log Kow for organic compounds [16,17]. As seen in Table 1, log Kfi^. for the training set showed a 
positive correlation with their log Kow (r = O.llA, p < 0.05). On the other hand, an inversely 
proportional trend was observed between log Kfi^- and water solubility (/- = 0.559, p < 0.05), 
indicating that the less soluble compound is more likely partitioned to SPME/PDMS than more 
soluble compounds. 

The molecular connectivity index conveys the degree of branching for the molecular structure [18]. 
The greater the degree of branching in a molecule, the lower will be the value of the molecular 
connectivity index. Table 1 shows that Kf\^. was well proportional to their molecular connectivity 
indexes with a correlation coefficient r = 0.743 (p < 0.05). Several researchers have indicated that the 
molecular connectivity indexes were significantly correlated to polarizability [12,18]. In this study, the 
molecular connectivity indexes of the training set were proportional to their polarizability (r = 0.676, 
p < 0.05). However, it is speculated that the branch molecule is less easily partitioned to SPME/PDMS. 

2.2 Empirical Models 

The best predictive model for Kfiv was obtained through the multiple regression analysis using a 
statistical software package SPSS 20.0 version (IBM, New York, NY, USA, 2011). Based on the 
stepwise regression with 95% confidence limits, the multiple regression analyses of was as follows 
(r = 0.99, ;?< 0.05): 



The above model explained 98% of the variance in the Kj^. data of the training set. Nirmalakhandan 
and Speece [12] developed a quantitative structure-activity relationship (QSAR) model for estimating 
Henry's constant. In their study, the best predictive equation consisted of the molecular connectivity 



log Z^,, = 1 .545 + 0.724O + 0.292 - 0.227 / + 0.001 ^5 



(4) 
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index, polarizability, and indicator variable of organic compounds. For the physicochemical 
descriptors of Equation (4), Table 1 indicates that the water solubility (r = 0.559) and the indicator 
variable (r = 0.583) were correlated to Kj^, with the lower correlation coefficients. If the water 
solubility and the indicator variable were excluded for Equation (4), the multiple regressions of 
Kfiv were as follows: 



log Z^,, = 1.447 + 0.7280 + 0.316 'j- 0.224/ 
log =1.502 + 0.9180 + 0.177 



(5) 
(6) 



For Equations (5) and (6), the correlation coefficients (r) were 0.989 and 0.984, respectively, which 
were similar to that of Equation (4). Figure 2 presents a comparison of Kfi^, for the training set with Kfiv 
predicted from Equations (4)-(6). Figure 2 shows that the agreement between the training set and the 
models predicted Kfi^. is very satisfactory. 

Figure 2. Comparison of Kfiv between the training set and predictive models. 
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Several researchers have used the linear solvation energy relationship (LSER) method to estimate 
Kfiv for SPME [19,20]. The LSER model consists of five descriptors of chemicals, i.e., molar 
refraction, polarizability, hydrogen bond acidity, hydrogen bond basicity, and McGowan characteristic 
volume (mainly representing London dispersion). As compared with the models developed herein, the 
physicochemical descriptors for the LSER method require a large data set. In particular, 
the physicochemical property data for the LSER method may be unavailable for the compound of 
interest [16,21]. 

To assess the validity of the above empirical models, Kf^. for the organic compounds of the test set 
was predicted using Equations (4)-(6). Figure 3 indicates a comparison of Kf\^, for the test set with Kf„ 
predicted using the empirical models. As shown in Figure 3, most of Kfi^, for the test set were greater 
than those predicted using the empirical models. The slopes of regression analysis shown in Figure 3 
implied that there was a systematic difference in the experimental measurement of Kf\^, between the 
training and test sets. For example. Table 3 indicates that the values of were 566 and 564 for 
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ethylbenzene and /?-xylene, respectively. In the study of Xia et al. [14], Kj\,, was 512 and 575 for 
ethylbenzene and /?-xylene which were not included in the test set. However, the correlation 
coefficients of K^,, were 0.813, 0.819 and 0.771 for Equations (4)-(6), respectively. For a practical 
purpose, the regression equation, consisting of the molecular connectivity index, polarizability, and 
indicator variable, can be an adequate model to predict Kj\v for SPME/PDMS. 

Figure 3, Comparison of Kj\v between the test set and predictive models. 
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Table 2. The training set of PDMS-water partition coefficients and 
physico-chemical descriptors. 



Compound 






O 


/ 






log Ko„ 


Benzene 


58 


2 


0.165 


1 


78 


1790 


2.13 


Toluene 


189 


2.411 


0.502 


1 


92 


526 


2.73 


o-Xylene 


485 


2.827 


0.839 


1 


106 


178 


3.12 


m-Xylene 


533 


2.821 


0.839 


1 


106 


161 


3.2 


jC-Xylene 


564 


2.821 


0.839 


1 


106 


162 


3.15 


Ethylbenzene 


566 


2.971 


0.839 


1 


106 


169 


3.15 


Cyclopentane 


712 


2.5 


0.733 


0 


70 


156 


3.01 


Methylcyclopentane 


1356 


2.894 


1.07 


0 


84 


42 


3.37 


Isopropylbenzene 


1412 


3.354 


1.176 


1 


120 


61.3 


3.66 


1 ,3,5-Trimethylbenzene 


1451 


3.232 


1.176 


1 


120 


48.2 


3.42 


1 -Methyl-4-Ethylbenzene 


1581 


3.382 


1.176 


1 


120 


94.9 


3.63 


Cyclohexane 


1592 


3 


1.07 


0 


84 


55 


3.44 


«-Propylbenzene 


1664 


3.471 


1.176 


1 


120 


52.2 


3.69 


1 ,2,4-Trimethylbenzene 


2183 


3.238 


1.76 


1 


120 


57 


3.63 


ter?-Butylbenzene 


2185 


3.661 


1.513 


1 


134 


29.5 


4.11 


2 ,3 -D imethy Ibutane 


2359 


2.643 


1.782 


0 


86 


22.5 


3.42 


2-Methylpentane 


3224 


2.77 


1.782 


0 


86 


14 


3.21 


3 -Methylpentane 


3270 


2.808 


1.782 


0 


86 


17.9 


3.6 
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Table 2, Cont. 



Compound 


Kf„ 


1 

1 


O 


/ 






log Ko^ 


1 -Methyl-3 -Isopropylbenzene 


3284 


3.765 


1.513 


1 


134 


42.5 


4.5 


trans- 1 ,2-Dimethy Icyclopentane 


3372 


3.207 


1.407 


0 


98 


33.9 


3.52 


1 -Methyl-3 -«-Propylbenzene 


3772 


3.882 


1.513 




134 


9.09 


4.67 


«-Butylbenzene 


3872 


3.971 


1.513 


1 


134 


11.8 


4.38 


1 ,2-Dimetliyl-4-Etliylbenzene 


3984 


3.8 


1.513 


1 






4.5 


i'ec-Butylbenzene 


4011 


3.892 


1.513 




134 


17.6 


4.57 


Isobutylbenzene 


4197 


3.827 


1.513 




134 


10.1 


4.68 


cw-l,3-Dimethylcyclopentane 


4289 


3.288 


1.407 


0 


98 


33.9 


3.52 


1 ,3 -Dimethyl-2-Ethylbenzene 


4345 


3.805 


1.513 


1 


134 


19.6 


4.28 


Methylcyclohexane 


4657 


3.394 


1.407 


0 


98 


14 


3.61 


cis-trans-cis- 1 ,2,4- 
















5621 


3.698 


1.744 


0 


112 


14.8 


3.94 


Trimethylcyclopentane 






trans- 1 ,2-Dimethylcyclohexane 


6638 


3.805 


1.744 


0 


112 


5.2 


4.01 


1 -Ethyl- 1 -methy Icyclopentane 


6831 


3.768 


1.744 


0 


122 


11.6 


4.05 


2,2-Dimethylpentane 


7349 


3.061 


2.119 


0 


100 


4.4 


3.67 


cis- 1 ,2-Dimethylcyclohexane 


7826 


3.805 


1.744 


0 


112 


6 


4.01 


«-Pentylbenzene 


8195 


AMI 


1.85 


1 


148 


3.37 


4.9 


2,4-Dimethylpentane 


8989 


3.126 


2.119 


0 


100 


5.5 


3.63 


2-Methylbutylbenzene 


9099 


4.365 


1.85 


1 


148 


12.7 


4.43 


2,2,3-Trimethylbutane 


9802 


3.944 


2.119 


0 


100 


28.9 


3.59 


2-Methylhexane 


10,202 


3.27 


2.119 


0 


100 


2.54 


3.71 


3 ,3-Dimethylpentane 


10,963 


3.121 


2.119 


0 


100 


5.92 


3.67 


3-Methylhexane 


11,146 


2.9 


2.119 


0 


100 


4.95 


3.71 


2,3-Dimethylpentane 


13,074 


3.181 


2.119 


0 


100 


5.25 


3.63 


1 ,2,4-Triethylbenzene 


16,253 


4.92 


2.187 


1 


162 


2.9 


5.11 


1 ,3,5-Triethylbenzene 


18,517 


4.914 


2.187 


1 


162 


0.014 


5.11 


2,2,3-Trimethylpentane 


21,205 


3.481 


2.456 


0 


114 


2.4 


4.09 


2,5-Dimethylhexane 


23,519 


3.626 


2.456 


0 


114 


9.2 


4.12 


2,2-Dimethylhexane 


24,504 


3.561 


2.456 


0 


114 


0.2 


4.16 


2-Methylheptane 


25,806 


3.77 


2.456 


0 


114 


7.97 


4.2 


4-Methylheptane 


27,274 


3.808 


2.456 


0 


114 


7.97 


4.2 


3-Ethylhexane 


28,370 


3.846 


2.456 


0 


114 


7.97 


4.2 


3-Methylheptane 


31,856 


3.808 


2.456 


0 


114 


0.792 


4.2 


2,3-Dimethylhexane 


33,749 


3.681 


2.456 


0 


114 


9.2 


4.12 


2,4-Dimethylhexane 


41,133 


3.664 


2.456 


0 


114 


9.87 


4.12 


2-Methyloctane 


45,267 


4.27 


2.793 


0 


128 


2.87 


4.69 


3 ,3 -Diethy Ipentane 


63,718 


4.243 


2.793 


0 


128 


2.9 


4.65 


3-Methyloctane 


66,682 


4.308 


2.793 


0 


128 


2.87 


4.69 


2,3-Dimethylheptane 


68,675 


4.181 


2.793 


0 


128 


3.37 


4.61 


3 ,3 -Dimethylheptane 


76,013 


4.121 


2.793 


0 


128 


2.9 


4.65 


3,5-Dimethylheptane 


78,829 


4.202 


2.793 


0 


128 


3.11 


4.61 


3,3-Dimethyloctane 


82,430 


4.621 


3.13 


0 


142 


2.9 


5.14 


2,5-Dimethylheptane 


84,142 


4.164 


2.793 


0 


128 


3.11 


4.61 


3 ,4-Dimethy Iheptane 


93,292 


4.219 


2.793 


0 


128 


3.11 


4.61 



Mw: g/mole; ^ WS: mg/L. 
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Table 3. The test set of PDMS-water partition coefficients. 



Compound 




X 


Q 


1 




Benzyl alcohol 


0.447 


2.580 


-0.323 




A A ^ / 1 a4 

4.29 X 10 


4-Fluorophenol 


0.525 


2.234 


-1.110 




1 ^ 1 a4 

1.25 X 10 


T\1 1 

Phenol 


0.661 


2.134 


-0.540 




o ^ o > , I a4 

8.28 X 10 


w-Cresol 


0.933 


2.545 


-0.323 




2.27 X 10 


Phenethyl alcohol 


1.318 


3.081 


0.014 




2 X 10 


3-Methylbenzyl alcohol 


1.479 


2.991 


0.014 




5 X 10"^ 


3-Chlorophenol 


2.042 


2.647 


-0.727 




2.6 X 10'* 


3 ,5-Dimethylphenol 


2.630 


2.956 


0.014 




4.88 X 10 


O T» 1 1 

3-Bromophenol 


2.884 


3.026 


-0.762 




/I 1 a4 

2.4 X 10 


4-Ethylphenol 


3.981 


3.106 


0.014 




AS X 10 


4-Chloroanume 


6.918 


2.676 


-0.022 




o A 1 a3 

3.9 X 10 


Phenyl acetate 


7.244 


3.023 


-1.430 




4.64 X 10^ 


Benzonitrile 


10.96 


2.384 


0.753 




2 X 10^ 


Acetophenone 


10.96 


2.865 


-0.605 




6.13 X 10^ 


4-Chloroacetophenone 


43.7 


3.342 


-0.672 




111 


Methyl benzoate 


A A ^ 

44.7 


2.977 


-1.430 




1 1 a3 

2.1 X 10 


Ethylbenzoate 


131.8 


3.565 


-1.093 




720 


A 1 * 1 

4-Chloroamsole 


T> /I A 

234.4 


3.036 


-0.390 




237 


Chloro benzene 


251.2 


2.477 


A A AO 

0.098 




/I AO 

498 


Bromobenzene 


323.6 


2.891 


0.063 




410 


T J 1 

lodobenzene 


537.1 


3.161 


0.692 




O /I A 

340 


Naphthalene 


676.1 


3.405 


-0.344 




31 


4-Chlorotoluene 


741.3 


3.095 


0.435 




106 


Propylbenzene 


1380 


3.471 


1.176 




52.2 


1 -Methylnaphthalene 


1819 


3.821 


-0.007 




25.8 


Biphenyl 


2344 


4.071 


-0.007 




6.9 




t 


WS: mg/L. 









3. Method 

3.1. PDMS-Water Partition Coefficients 

The training set of PDMS-water partition coefficients for alkanes and aromatic hydrocarbons was 
compiled from the literature [13] published by Dr. Pawliszjoi who invented the SPME technique. As 
shown in Table 2, the values of Kfw (n = 61) exhibit a wide range from 58 (for benzene) to 82,430 
(for 3,3-dimethyloctane). In addition, the organic compounds of the training set have a wide range of 
water solubility, e.g., 1.4 x 10 mg/L for 1,3,5-triethylbenzene, and 1.79 x 10 mg/L for benzene. 
The test set of PDMS-water partition coefficients (n = 26), not included in the training set, was 
retrieved from the study of Xia et al. [14]. Table 3 indicates that the magnitude of Kj\v values for the 
test set, consisting of alcohols, benzenes and phenols, were between 10"^ and 10^. Kj^ for the training 
and test sets were experimentally obtained using the SPME/PDMS devices at 25 °C. 
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3.2. Physicochemical Descriptors 

In this study, the polarizabihty and molecular connectivity index were used to represent the 
interactions between the analj^e and PDMS coating. As shown in Equation (7), the polarizability 0 
was determined by the addition of atom/bond contribution factors [12,22]. 

0 = ^F.x(number of atom /bond) (7) 
where is the contribution factors indicated in Table 4. 



Table 4. Contribution factors for the calculation of polarizability. 



Atom/bond 


Contribution factor 


Atom/bond 


Contribution factor 


Carbon 


0.577 


Iodine 


0.407 


Hydrogen * 


-0.120 


Fluorine 


-0.570 


Oxygen 


-0.825 


Cycle 


-0.952 


Hydroxyl 


-3.701 


Double bond 


-0.859 


Chlorine 


-0.187 


Triple bond 


-0.109 


Bromine 


-0.222 







Data was obtained fromNirmalakhandan and Speece [12]; * Attached to carbon atoms only. 



The molecular connectivity index is a topological descriptor of molecular structure based on a count 
of skeletal atom groupings of a chemical compound [18,23]. For each atom of a molecule, a 5 value is 
assigned as the difference between the number of valence electrons and the number of hydrogen 
atoms attached to that atom. In this study, the first-order connectivity index, ^x? was calculated as 
follows [9,11,12]: 

'z-tmr;-' (8) 

where 6, and 6, are the 6 values of two adjacent atoms i and j, respectively; and n is the number of 
bonds in the molecule. 

An indicator variable, /, was used to differentiate between compounds on the basis of their ability to 
take place in hydrogen bonding [24]. / was assigned a value of 1 for all compounds containing an 
electronegative element, such as oxygen, nitrogen, and halogen etc., attached directly to a carbon atom 
holding a hydrogen atom. In addition, the value of / was assigned 1 for acetylinic compounds and 
aromatic compounds with partially substituted hydrogen atoms. For the other compounds, / was set 
equal to zero [12]. 

4. Conclusions 

From the single-parameter model analyses of the training set, the PDMS-water partition coefficients 
were proportional to the molecular connectivity index, polarizability, molecular weight, and 
octanol-water partition coefficient of organic compounds, while an inversely proportional trend was 
observed between Kjw and the water solubility as well as the indicator variable. Kjw was significantly 
dependent on the polarizability of the organic compounds. Based on the results of multiple regression 
analyses, several correlations for Kj^ were developed using the polarizability, molecular connectivity 
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index, indicator variable and water solubility, indicating a good agreement {r > 0.11 \) between the 
training set and the predicted Kj^- In addition, the correlations were able to adequately predict Kj^^ for a 
test set which were not included in the correlation development. The empirical model developed in this 
work is more versatile than current available correlations. In order to broaden the predictive abilities, 
however, future work should be conducted to calculate Kj^ using the empirical model for a variety of 
organic compounds. 
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